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TECHNICAL NOTE 3176 


WALL INTEIRFERENCE IN WIND TUNNELS WITH SLOTTED AND 

POROUS BOUNDARIES AT SUBSONIC SPEEDS 

By Barrett S. Baldwin, Jr., John B. Turner, 
and Earl D. Khechtel 


SUMMARY 


Linearized compressihle-flow analysis is applied to the study of 
wind-tunnel-wall interference for subsonic flow in either two- 
dimensional or circular test sections having slotted or porous walls. 
Expressions are developed for evaluating blockage lift interference. 


INTRODUCTION 


In solid-wall wind tunnels the effects of blockage severely limit 
model sizes that can be tested at hi^ subsonic speeds; in fact, the 
model must become vanishingly small as sonic speed is approached. It 
has been demonstrated that if the walls are ventilated (e.g., slotted 
or porous) then blockage is reduced and much larger models can be 
tested. However, wall -interference effects, although reduced, still 
exist and must be evaluated in order to correct the wind-tunnel data to 
free-air conditions. 

It is the objective of the present investigation to analyze two of 
the principal wall- interference effects, blockage and lift Interference, 
for two- and three-dimensional subsonic flows in ventilated test 
sections, where blockage refers to the mean Incremental velocity induced 
in the vicinity of the model by wall Interference and lift interference 
is the mean upwash so induced. In the three-dimensional case it is 
convenient to perform the analysis for a circular test section. The 
results obtained for the circular test section may be applied to a square 
test section of equal cross-sectional area since the wall interference 
at the center of the tunnel should be relatively insensitive to such a 
change in the shape of the wall. 
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SYMBOLS 

factor In Fourier integral transform of (p* 
slot width of slotted wall (see fig. l) 
wing span of model wing 

constant factor in nonlinear term of boimdary equation 
Fourier integral transform with respect to x of cp 
dummy variable of Fourier transform 
half tiinnel height 

modified Bessel function of the first kind and order zero 

modified Bessel function of the first kind and order one 

modified Bessel function of the second kind and order zero 

modified Bessel function of the second kind and order one 

slot constant, — — Zn 

slot separation of slotted wall (see fig. l) 
lift on the model 
free-stream Mach number 

T IT 

parameter proportional to size of two-dimensional model, 

parameter proportional to size of _ three-dimensional model, TjiU 

coordinate in the direction of the outward normal to the wall 

dummy variable of Integration 

porosity parameter 

cylindrical coordinates 

free-stream velocity 
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u,v,w 


U*J V*,-W* 


^ /V ru 

u,v,v 


Aii,Av,Aw 


W 

X 

x,y,z 

a 

P 

r 

■3> 


<p 

9i 

cp* 

69* 


p 


Te 

Tr 

I 


pertTirbation velocity components In the x, y, z directions, 
respectively 

additional velocity components due to the presence of the 
walls 

additional velocity components having rapid spacewise 
variation near the walls 

additional velocity components at the position of the model 
due to the walls 

complex velocity in the y, z plane 

complex variable equal to z + iy (physical plane) 

Cartesian coordinates 

dummy constant in limiting process 

circulation 

complex velocity potential 

total perturbation velocity potential, 9i + 9* 

approximate perturbation potential due to model in free air 

additional perturbation potential due to tvinnel walls 

additional wall-interference potential arising from non- 
linear term in boundary equation 

free-stream density 

cross-sectional area of two-dimensional model 
volume of three-dimensional model 

function of X, equsl. to C + iq (transformed plane) 
Cartesian coordinates in tremsformed plane 
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ANALYSIS 1 

General Statement of the Problem 


The effect of-the tunnel walls on the flow around a model, in the 
case of ventilated walls, can be calculated using the same basic method 
as that used in reference 1 for the closed -wall case. As in reference 1 
the analysis is based on the linearized equation of subsonic compressi- 
ble flow 




^ 

Sy^ 


( 1 ). 


where cp is the pertxrrbation velocity potential of the flow in the 
tunnel. 

Let qp = cpi + Cf*-, where Cfi is the potential of the flow about the 
model in free air and 9 * is the potential of the additional flow due 
to the presence of the walls. 

If cpi is taken to be a known solution of equation (l) which 
approximates the true free-air potential at points far from the model, 
9* can be calculated from the fact that the sum 9i + satisfies a 
known boundary condition at the wall. Since the values of 9i at the 
wall only are used, any inaccuracy in the value of 9i near the model 
should not affect the calculation of- - 9* appreciably. 

The primary objective in this procedure is to estimate the change 
in stream conditions caused by the walls at the position of the model. 

It is assimied that the velocity conqjonents derived from 9* 8-re con- 
stants near the model which can be subtracted from the stream velocity 
to obtain the equivalent free air stream velocity. Thus, 

Au = at x=y^z=0 

ox 

is the blockage correction, and 

Aw = at x=y=z=0 

oz 

is the upwash correction in the three-dimensional case. 


Boundary Conditions 


In this section a single expression approximately representing the 
boundary conditions of solid, porous, and slo^tted walls and an open jet 
will be developed. 
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Let X "be the coordinate in the direction of the free stream 
and n the coordinate in a direction perpendicular to^the x . direction. 
Consider a wall which is perpendicular to the n direction (i.e., 
parallel to the free stream) . If the wall is solid, the condition of 
no flow through the wall can he expressed as 

— = 0 at the wall 

on 

In the case of an open Jet there is no pressure drop across the jet 
boundary so that there is zero perturbation pressure at the boundary. 
With a disturbance in the stream this boundary does not remain paraiJLel 
to the free stream. However, for convenience, the condition of zero 
perturbation pressure is imposed at a surface parallel to the free 
stream and coinciding with the jet boundary far upstream of the 
disturbance (see ref. 2). Also, for convenience, this surface can be 
called an open wall and the bounda 2 ry condition can be expressed as 

— = 0 at the wall 

Sx 

In reference 3 an average boundary condition for a porous wall is 
derived. The average velocity normal to the wall is assumed to be 
proportional to the pressure drop through the wall, a linearized 
approximation to viscous flow throusfh a porous medium, and the pressure 
outside the wall is assumed equal to the free-stream pressure. This 
leads to the boundary equation 

+ i ^ = 0 at the wall (2) 

5x R Sn 


The quantity R is a porosity parameter defined by 

^ =£U ^ 

R Sn 


(3) 


where 

Ap pressure drop through the wall 

P stream density 

U stream velocity 

The quantity PU/R can be determined, experimentally by measuring the 
mass flow and pressure drop throu^ a sample of the wall lander con- 
ditions corresponding to zero stream velocity. 

Porous walls to which eqToations (2) and (3) axe applicable will 
henceforth be referred to in this report as ideal porous walls. 
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An approxi ma te boundary equation for a slotted wall is derived in 
Appendix A. The pressure at the slo.ts is assumed constant and equal to 
the free-stream pressure . The resulting uniform boundary condition is 


+ K — = 0 at the wall 

ox Sx^n 

where K is related to the slot geometry by 



(h) 

(5) 


Slotted walls to which equation (U) is applicable will henceforth be 
referred to in this report as ideal slotted walls. 


Solutions for wall interference based on equations (2) and (4) can 
be obtained in one calculation by combining them in the form 

— + K ^ = 0 everywhere at the wall (6) 

9x 3x9n R Sn _ . 

Thus, wall-interference solutions based on equation (6) contain, as 
special cases, those of the closed wall (K -> » or 1/R -5> »), ideal 
porous wall (K = 0), ideal slotted wall (i/R =0), and open jet 
(K = 0 and 1/R = 0) . Furthennore, equation (6) can be assumed to 
describe a slotted wall having mixed potential and viscous flows in the 
slots. In that case the porosity parameter R can be determined 
experimentally, as it is in the case of a porous wall, by measuring the 
mass flow . for a given pressTire drop through a_ sample of the wall. 


If it is found that a nonlinear relationship between pressure drop 
and mass flow exists, it may be necessary to add a term of the form 
f(S9/Sn) to equation (6) . This case is discussed in Appendix B. 


In addition to the foregoing interpretations of equation (6) , an 
interpretation identifying it with slotted walls with tapered slots 
(l and a functions of x) and potential flow in the slots is possible. 
This case is discussed in Appendix C where it is found that Instead of 
representing viscous effects, the parameter R is related to the taper 
i3y 


where K(x) is the same as in eqiiation (5)« 


1 

R 


dx 


i X = 0 


(7) 


Blockage in a Two -Dimensional -Flow Tunnel 


Under the assumption of infinitesimal model size, the blockage cor- 
rection will be calculated using equation (6) as the boundary equation. 
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Let X be the coordinate in the free -stream direction and y the 
coordinate in the direction perpendicular to the walls. Let h be the 
half tunnel height .so that the walls are at y = -h at y = +h. 

In these coordinates eq,uation (l) becomes 


and equation (6) 




= 0 


becomes 

3<p 


± K 


ojray 


1 a<p'\ 


y = ±h 


=0 


(8a) 

(8b) 


The t signs on the second and third tenns are required because at 
the upper wall, n = +y and at the lower wall, n = -y, n being the 
coordinate in the direction of the outward normal to the wall. 


If 9 is replaced by 


9i + 9*, 






equations (8a) and (8b) yield 

= 0 (9a) 


0 


^ ± K ± i ^ 

Sx ^x8y E 8y /y _ 


Sx 


± K 


± i (9b) 

dxoy R dy/ y = ±h 


These two equations are sufficient to determine 9* when 9^, the dis- 
turbance due to the model in free air, is known. 

As in reference 1, the disturbance due to the model at zero angle 
of attack in free air is approximated by a two-dimensional doublet 
which can be expressed as 

- ■ s (^) <« 

The reasoning behind the choice of the doublet is as follows: The 

source-sink distribution representing a nonlifting model contains the 
same total sink strength as total source strength, so that the distant 
flow field would not resemble that of a single source or sink. The 
center of gravity of the source distribution would lie forward on the 
model compared to the center of the sink distribution which would be aft. 
Hence, there would be a dipole or doublet effect of the model at a 
distance. Other higher-order effects would be much less than that of 
the doublet at large distances. The constant m^ is related to the 
size of the model with sufficient accuracy for present purposes by 


V 
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vhere Tg 


me = 


TeU 

P 


is the area of the model in the 


xy plane . 


( 11 ) 


It is convenient to solve equations ( 9 ) by the Fourier transform 
method. The Fourier transform of is defined by the relations 

00 

Gi(g,y) cpi(x,y)e^e^ dx (l2a) 

-00 





Gi(g,y)e"^^ dg 


(12b) 


Analogous expressions defining the Fourier transforms of 


q> 


are 


00 


G*(g,y) = i 

2it , 

/ 1>*(x,y)e^e* dx 
-00 

(13a) 

00 

^ (x,y) = y 

-00 

0 dg 

(13b) 


With the substitution of equation (l3b) into (9a) and with an interchange 
in the order of differentiation and integration there results 

00 00 

f ^ G*(g,y)e‘^®^ dg + J ^ G*(g,y)e"^sx ^ q 


Performing the indicated x differentiation and collecting the two 
terms under the common integral yields 



This integral will be zero if 






= 0 


(l4a) 


Substituting equations (l2b) and (l3b) fnto equation (9b) and 
proceeding as before yields 
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( . n* ~ OG + 1 SG A f , „ . „ . 1 5Gi ^ 

\ Sy R ^/y ^ V ^ Sy R oy /y = ±i^ 

(.ihh) 

Eqmtions (l4a) ajid (l4b), involving the transforms Gi and G* of 
*?i and cp*> can now he solved in the place of solving equations (9a) and 
(9b) . The general solution of equation (l4a) is 

G*" = A(g) cosh (Pgy) + E(g) sinh (Pgy) 

Substituting the second term of this expression into equation (l3b) 
yields a term of cp* having an odd y dependence. From symmetry 
considerations it can be seen that cp* should have even y dependence. 
Therefore, B(g) must be zero and, hence, 

= A(g) cosh (pgy) ( 15 ) 

The unknown factor A(g) will be determined by substituting equation ( 15 ) 
into equation (lltb), but Gi must be known for this purpose. 

In order to find Gi, equation (lO) is substituted into equa- 
tion (l2a) and there results 


(s,v) = g 


From reference I 4 -, it is found that 

00 


i r L 

2jt J f 


dx = e' 

2P ly| 


Iplyl 


GO 

^In reference It, G(g) is defined as G(g) = j F(f)e^^”^S ^f, in 

-CO 

the tables of F(f) versus G(g), F is usually given as a function 
of p where p = 12jtf. If P(p) is the function of p given in the 
table, it follows that F(f) = P(p) = P(i2rtf). With the substitution. 


X = 2«f, it is found that G(g) = 




P(lx)e'*^ dx. Thus, G(g) in 


the tables is the Fourier transform, as defined in equations (l2) of 
the function of x which resiolts from replacing p by ix in P(.p) . 
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vhere the syuibol | J denotes absolute value. Differentiating both 
sides of this with respect to g yields 




so that 


Gi = ^ i e'lslPkl 

^ (g| 

Substituting equations (15) and (17) into (lijb) yields 


(17) 


A(g) 


- ig cosh (Pgy) T Ig K Pg slnh (pgy) ± ^ Pg slnh (Pgy) 


y = ±h 


8 

^ ^ Tli 


ig e 


- g( P[y 


Tig K 




glP|y 



Making the indicated substitutions (y = ±h) and solving for A(g) yields 
the single complex equation ... 




(lf|j--lgK6 +£)e-|e|P» 


cosh (pgh) + Kpg sinh (p^) + i | sinh (Pgh) ] 


J 


(18) 


Substituting eqviatlons (15) and (18) into equation (l3b) yields 




( 


1 -T^ — ig Kp + £ ^e cosh (Pgy) 

R , 


cosh (pgh) + KPg slnh (Pgh) + i slnh (Pgh) 

R 


<3-g 
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Updn separation into real and Imaginary parts and making use of the fact 
that an Integral from -00 to +« is twice the value of the integral 
from 0 to 00 if the integrand is an even function and zero if it is an 
odd function, this becomes 



With the substitution, q. = Pgh, this becomes 




t 
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The integrals indicated in equation (19) have not been found in 
closed form. However, if the derivatives of cp* are desired at a 
small number of points for a small number of values of K/h and 3 /r, 
the derivatives can be taken in the integrand and the resulting inte- 
grals evaluated numerically or graphically. 


The quantity of primary interest is Au, the value of Sq>^9x 
evaluated at x = y = 0. Performing the indicated differentiation on 
equation (19) and setting x = y = 0, the expression 



is obtained for the additional stream velocity at the position of 
the model due to the walls . 


Solid wall .- Letting K — > » or 1/R — > » in equation (20) yields 


m. 


Au 


/- 


e~'3q^ 


jt nie 

dq = — 


27tP^h^J sinh (q) 24 ^®h^ 

h ° 

Ideal porous wall .- At K = 0 equation _ (20) becomes 

“e 


(21) 


Au 


h = 0 


^ '^o cosh (q) 


/ \2 

" cosh (q) - f sinh q 
J- r i_ / _ \ ”i r/' b\ . , /■ 


+ If £ ) sinh (q) 


e"\dq (22) 


t 


ft 


Ideal slotted wall .- With 1 /r = 0 we have 

- f q) e"^\dq 


Au 


P = o 


me 

jtp^h^ 


/ 




Open jet .- For both K = 0 and 1/R = 0 


Au 


me 


i-o 

I- 




J 1 + e"' 


^dq 

•aq 


n mg 


(23) 


(24) 


» 
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In the general case (eq.. 20), the value of Au lies between the 
values for the solid wall and the open jet, and Zfli can be made arbi- 
* trarily small by an appropriate choice of K/h and p/R. Figure 2 is a 

plot of the values of K/h and P/r at- which Au = 0. These values 
were computed numerically near the ends of the c^lrve and interpolated 
in the middle. Eq.uations (21), (22), and {2k) are in agreement with 
the results of reference 3* Figure 3 is a plot of equation ( 23 ), show- 
ing the variation of blockage factor with slot parameter for the two- 
dimensional-flow, ideal, slotted tunnel . 

Since the effect on the blockage correction of letting P approach 
zero is the same as letting 1/R approach zero, it can be concluded 
that as the stream Mach number approaches unity, the blockage correction 
factor for an ideal porous wall approaches that of an open jet. Simi- 
larly any effect of viscosity or taper of a slotted wall described by 
equation (6) would be suppressed at near sonic speed so that the block- 
age correction factor would approach that of an ideal slotted wall. 


Blockage in a Circular Tunnel 


Again, the blockage due to the wall interference of a very small 
model will be calculated using equation (6) as the boundary equation. 


Let X be the coordinate in the free-stream direction and r the 
cylindrical coordinate perpendicular to the x direction. Let rQ tie 
the tunnel radius so that the wall is at r = Tq. Using these coordi- 
nates, equation (l) becomes, in the case of rotational Invariance, 




^ 1 S 

Sx^ r 3r 



= 0 


( 25 ) 


and equation (6) beccxnes 


( 


^ + 
Sx 


K 


. 1 


Sxdr 


Sr / 


= 0 


r = ro 


( 26 ) 


Then 9* must satisfy the equations 

^ ^ + i A fr ^') = 0 (27a) 

Sx^ r Sr V Sr / 


(is! , K * I , K * i (aTb) 

V Sx SxSr R Sr / r = ro \ 5x SxSr R Sr ^ 

(See fig. 1.) 
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Again, from reference 1 the free-alr solution is approximated by 

"" (28) 


'Pi = 


W (x2 + p2^)3/S 

The constant mp is related to the size of the model by 


m^ = U Tj. 


(29) 


where is the volume of the model. Correction factors which take 

into account the shape of the model and the presence of shock waves are 
discussed in reference 5. 

Substituting, as before, the Fourier integral expressions for 
and 9 i into equations (27) yields 


(- 


ig G - K ig 


Sr R 8r y 


(30a) 


5 Gi 


-ig Gi - K ig ^ + 
or 


1 SgA 

R SrZ 


k* = r. 


(30b) 


The general solution of the first of equations (30) is 
G* ■= A(g) Iq (Pgr) + B(g) Ko (Pgr) 

where Iq and Kq are the modified Bessel functions of zero order of the 
first and second kinds, respectively (see ref. 6). Since Kq (Pgr) has 
a singularity at r = 0 which would lead to a singularity in q>* at r = 0, 
B(g) must be zero and, hence. 


G* = A(g) Iq (Pgr) 


(31) 


Substituting equation ( 28 ) into the Fourier Integral exqiression of 
9 i yields 

00 

mr 1 r X 


Gi = ^ 


Iht _ J (x2 + p2r2)®/2 

From reference k it is found that 


igx 

e dx 


(32) 
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According to eq,uations (12), the inverse relation 


follows. Taking the derivative with respect to x of both sides of 
this equation yields 


00 

J' (-Ig) ^ Kq (Pr[g|) 


e-^Sx ^ ^ . 


(x^ + 

From this, according to equations (l2), 

00 

— f 7 -^ ^ ■ „ 3 7 2 - dx = ig 1 Ko (3r|g| ) 

2k J (x^ + p^r^) / 


so that 


A(g) 


Gi = ^ Ig Ko (Prlgl) (33) 

Substituting eqizations (31) and (33) into equation (30b) yields 
-ig Iq (Pgro) -ig KPg Ii (PgTo) + ^ Pg Ii (Pgro)J = 


m. 


4k2 


ig 


-ig Kq (Prolgl) -ig K(-Plgj) Kj_ (Prolgl) + 


I (-3lgl) Ki (Prolgj) 


the derivatives of Iq and Kq being found in reference 6 . Upon 
solving for A(g) , substituting the result into equation (l3)# and 
substituting q for prQg there results 



cp = 


ON 


lOl’M 

€0 

a ? r 

r. — 1 . 

Ki(q)lo(q) + Ko(q)li(q) 

loCgys ('g-')a aa 

4i(P^ro^ 

R It 

Io(<l) + ^ <lli(<l) 
L ^0 

]' ^ [5 


IF 


Ki(q.)lo(q)+Ko(4)lo(q)+ ^ Ko(q)li(q)- ;^')Ki(q)li(q)] Iq sin 


\HAO/ 


T 


K T.^„^T 

-uj_i •' ■ • 


.. T. f r . 

J 


£ 


T, 

•**JL.\ U/ 


Differentiating vith respect to x and setting x = r »* 0 yleMs 


(31^) 


M = 


utj* 


■»£d3-wi 3 


’ iKL(a^I^(a) + a 

-Ki(q)lo(q) + Ko(q.)7i(q.) 

- 1 

fa^ rrC- + ~'^Kn f q)It (a)l 
V -^07 -p^y ~'~ J 

[■ Q^dq 


K ■ . 1 

I«f + a — Ii (a^ 

- Tq *--j 

2 

+ 

£ Ii (a) 
K * • - 

P 


Closed va]-l.- Setting 1/R — >« in equation ( 35 ) yields 


( 35 ) 


An 


= 




p J-iU) 2n-r/p® 

R 


( 36 ) 


Ideal porous vaH.- At K = 0 equation ( 35 ) becomes 
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Au 





€0 

1 / 


' K^(<l)lo(9) - ^ Ki(q)li(q) 

q^dq 


Io(<l) 

2 

+ 

^ Ii(cl)1 
Lr J 

2 


(37) 


Ideal slotted vail .- With (i/r) = 0 we have 


Au 


mr 2 

00 

r 

Ko((l) - CL — Ki(q) 

L I'o J 

q^dq 

4«3®rQ® 

3 ^ 0 c 

# = 0 

) 

Io(l) + <1 JL Ii(q) 

L ro J 



Open Jet .- For both (i/r) = 0 and K = 0 



00 

mr 2 r Kp(q) ^ 0.63 mr 


(38) 


(39) 


The values of Au in equation (35) lie between the values for the 
closed wall and the open jet. Figure 4 is a plot of values of K/ro 
versus P/R at which Au = Oj the shape of the curve was calculated near 
the ends and interpolated in the middle. A graph of equation ( 38 ) appears 
in figure 5 showing the variation of blockage factor with slot parameter 
for the cylindrical, ideal slotted tunnel. 

Again, letting 3 0 has the same effect on the blockage cor- 

rection as letting (i/r)- 5>0, so that near sonic speed the ideal porous 
wall should act like an open jet and the slotted wall should act like 
an ideal slotted wall. 


Lift Interference in a Circular Tunnel 


The upwash correction will be calculated using the infinitesimal 
model size approximation and the approximate boundary condition of 
equation (6). 

Let X be the coordinate in the free-stream direction, z the 
coordinate in the direction of lift on the model, and y the remaining 
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Cartesian coordinate. The cylindrical coordinates r ajid 0 are related 

to the Cartesian coordinates hy r = + and z = r sin 0. It can 

be seen from reference 1 that the appropriate free air solution Is 


«Pi = 


h-fl 



X ■ sin 0 


(40) 


which Is the potential of a horseshoe vortex having Infinitesimal span. 
The fact that the actual span Is finite introduces higher-order terms 
which are negligible at distances large compared to the size of the 
model. Here Tb is related to the lift on the model by 


L = pUFb 


(41) 


The Fourier transform with respect' to x of cpi cannot be found. 

An arbitrary parameter a will be introduced into the potential so that 
the Fourier transfonn. of a related function can be. found. Let this 
function be 


£b + P^r^ 1 -oaJ sin 0 

4ir \ a 8x / r 


(42) 


so that 


lim ’ = ^Pi 
a, — > o 


Then a will be eliminated from the resulting cp’*’* by taking the limit 
with a at zero . 

From reference 4 it is found that the Foiirler transform of 


_ + o2.^ cxgr Ki (Pr</g^ + a?) 

e P ^ is ■ — . By the use of an Inversion, 

It*/ g^ + 

differentiation, and reinverslon, as in the derivation of equation (33)^ 


It is found that the transform of — + 3 r^ 

dx 


is 
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igaPrKi ( 3r ,/g®+ a^) 


so that the transform of tPx' is 


(a+ig) 


sin 6 


(i*3) 


From physical reasoning, it can he seen that the sin 6 factor 
of cpj^ will result in the same angle dependence for cp*, so that sin 6 
will be a factor common to cpi * , 9*' and their Fourier transforms Gi ' 
and G* . With cylindrical coordinates, eq.uation (l) is 


r drV 3r/ r® 30^ 


In the case of sin 0 dependence, this becomes 


r Sr \ 5r / 


and the boimdary condition remains the same as in equation ( 26 ) . Sub- 
stituting the Fourier integral expressions for 9i* and 9* in these 
equations, as before, yields 


r dr Sr 


■) 


_ G*' =0 


(1^4a) 


(- 


igG** - Kig 1 


Sr R Sr A’=- 


ML ) = - (^-igGi ’ -Kig Mil -h 1 Mil) 

Sr ^r=r ^ Sr R Sr ''x=r 


(44b) 

The solution of equation (44a) vhich has no singularity at r = 0 is 


G* = A(g) Ii(3rg) sin 0 


(45) 



tdINxi 


Substitrutlng this and eq^uatlon (^3) into equation (4it-b) 3 rleldB 


ro 

o 


A(g)| - iK + iK K jl^.g^PS) _ ie K i Il(P£ Qg l + M T.^Rt- 1 nir, fl - 

■ ~ “ ■ ~U'' — tJO' j — " ■' 


3”0 


R To 


^ yr-KxP roVg2+a^ R 


TC-( n-K- 

■“U\I-*U'V O / J DJ-u. 




iui- cjuuB uj. ouu±ng ±ii wae iategraj. expresBion lor ^ ana i}aKing Tine ximix 
of q>*-^ ae a -» 0 yields 


r 


VM 

/ 


Fb 


Ihtrr 


sin 0 


__ P 2 n 

R nj 


Ki(q)lo(q) + Ko(q)li(q) 

_ 

I-I 

H 

ain 0 cos i 

^3roy 

1 da 

r * 


6- ■' — Vi(q) + <l^~ \ I.(q) 
LV ^0/ Vfo/ 

J W 

ii(q) 

/ ^^2 
- al-fa^ 

■*• U' 


2 

tt 


If V ^ . . 

— 

1 r / \ 1 / \^ 

■'i^Mi)JLV^-^i('i) +aj 5 lc(q) J+(|) 

1 r 1 

Ki(q)+qKb(q) ii(q)-qio(q)J 



^ ^ (IJ 

ii(q) - q io(q) 

r 




i j ® sin I ^ J dq 

\ ^ \ r--0/ 






f 
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n 


i n 


Jtteplacing ij. (qr/ro) tiy ItB power Berlee expansion, differentiating with respect to z, 
which is equal to r sin 6, and setting r = 0 yields 



00 

1_ £ 1 /" 

[Ki(q) Io(q) + Ko(q) Ii(q)] qS cot 

'(ft)- 

4rtro^ 

r=o 

R n^o qS 

Uiii. 


+ 

?i[ 

ii(q) - qio(q.) 

2 Jt 


/ 


2 ' 

Tn’l - n fl _ 2L ^T. E. T fn rw. ('n'l ^ T. T 


1 ~ — '\lifa') + a — Ir,(a) T + 
V ^o/ -'*■ * To j 


Infa^ - Q I^fa^ 


1 C^)'^ 

for the upwash at r = 0. 


(47) 


When X " 0 we have 


Pb 

00 

1 f 

[Ki(q) Io(q) + Ko(l) ^i(Q-) 1 ci®dq 



4rtro^ 

R J( ^ ^2 


ii(<i) - a io(<i) 

r. 


(M) 


for the value of upwash at the position of the model due to the walls. 


Solid wall.- Letting 1 /r m yields 


i 

R 


r-b_ 


00 


P 
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Ideal porous vail.- At K = 0 equation (48) becomes 


Awl 


Fb 


K 


4jrro^ 


■ =0 


1 _ P i 

R jt 


[Ki(q)l„(q) + Ko(q)li(q) ] q^dq 


( 49 ) 


q.^Ii^(<l) + 


f § 1 

Vr/ 


[Ii(q) - qIo(qL)l‘ 


Ideal 'slotted, vail . - When .1/R is set equal to zero in equa- 
tion (48), a limiting process is required to .obtain the correct result 
at (P/r) = 0. The result of this process is 


p 

R = ° 



Open jet. - Letting K = 0 


in equation ( 50 ) yields 


( 50 ) 



Fb 

4«^o^ 


The value of Aw in the general case of equation (48) lies between 
the values for the. closed wall and the open jet. 


CONCLUDING REMARKS 


A method of evaluating wall Interference of partly open walls involv- 
ing mixed potential and viscous flows has been presented. Expressions for 
blockage and lift interference for both slotted and porous walls have been 
derived. Some new details of the method may prove useful in other theo- 
retical treatments of this type of problem. 

The results of the analysis indicate that near sonic speed, the 
blockage correction for an ideal porous wall approaches that of an open 
jet. Similarly, any linear viscous or taper_effect of a slotted wall is 
suppressed near sonic speed, so that the blockage correction approaches 
that of an ideal slotted wall. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Eield, Calif., May 29, 1953 
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APPENDIX A 

DERIVATION OF THE BOJNDARY EQUATION 


In this appendix an approximate smoothed or average boundary equa- 
tion for a slotted wall will he derived. 


An ideal slotted wall has zero perturbation pressure at the slots 
and zero normal flow at the strips. These conditions can be expressed 
as 


^ = 0 at the slots 

Sx 


^ =0 at the strips 

Sn 


(Al) 


When the slot spacing and model dimensions are small compared to 
tionnel dimensions, the perturbation flow can be separated into a rapidly 
vai^ng and a relatively mlform part so that the two parts can be 
investigated separately. It will be shown that the effect of the rapidly 
varying pant can be replaced by a condition on the relatively imiform 
part. 


Let 9 , u, V, and w represent the rapidly varying part of the flow 
field and cp, u, v, and w the remaining part of the perturbation flow. 
For a plane wall at z = h equations (a 1) require that 


u + u = 0 at the slots 

w + w = 0 at the strips 


(A2) 


In addition 


u = V = w = 0 far from the wall 


(A3) 


In order to solve for 9 , use can be made of the fact that u, v, 
and w are nearly constant at the wall compared to u, v, and w, so 
that u and w can be considered constant in equations (A2). 

Since the slots lie along the x direction, u is nearly constant 
in the x direction, so that Bu/Sx can be neglected compared to 

and Sw/3z. As in slender airplane theory, this leads to a two- 
dimensional crossflow for which 
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(A4) 


Then $ can be of the form Q(x) f(y, z) where Q(x) is a slowly- 
varying function of X, and f(y, z) satisfies the two-dimensional 
Laplace eq.uation. Then, 


u 



f(y,z) 


and since u is equal to -u, a constant, at the slots, f(y, z) must be 
constant at the slots or 


V 




at the slots 


This equation can be used to replace the first of equations (A2), 
and altogether there results 


V = 0 at the slots 
w = -w = constant at the strips 
^ s= V = w = 0 far from the wall 


(A5) 


assuming that 5 q/Sx Is not absolutely zero. 

Equations (aU) and (A5) can be solved using the conformal trans- 
formation technique. 


Let the wall at z = h be slotted periodically from y = -00 to 
y = +00 with the center of a slot at y = 0. Let a be the slot width 
and 2 the slot separation (see fig. l) . It is sufficient to consider 
only one period of the periodic flow configuration so that solid boiinda- 
ries can be placed at y = ±(z/ 2 ) and attention confined to the region 
between them. 

Let X = z + iy be the complex physical plane and | = C + in the 
transformed plane. Let 9 be the canplex velocity potential such that 
W » w - iv, the complex velocity, is equal to d$/dX, and $ = real part 
Then eq-uation (a4) Is satisfied if 9 is any analytic function of X. 
The boundary conditions (A5) are satisfied by finding the analytic func- 
tion I (X) which transforms the boundaries in the X plane into a con- 
figuration for which the flow field with the desired flow at the bound- 
aries can be foimd. 
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In the X plane the region under consideration is from z = -» 
to z = h and from y = -(z/ 2 ) to y = +(z/2). The slot lies on the 
line X = h + iy from y = - (a/2) to +(a/ 2 ) . (See fig. 6.) The 
two half strips lie on the same line from y = -(l/2) to y = -(a/2) 
and from y = (a/ 2 ) to y = (1/2). The remainder of the solid boundary 
lies along the line X = z — l(z/ 2 ) from z = -oo to z=h and along 
the line X = z + i(z/2) from z = -« to z = h. 

The transformation which will place the origin at the wall, and the 
domain \mder consideration in the ri^t half -plane is 


61 = h - X (a6) 

The effect of the transformation on the positions of the boundaries is 
obtained by substituting the equations representing the boundaries into 
equation (a 6) • Thus, the line X ■= h + iy becomes + iTii = -iy 

in the li plane and the slot lies on this line from qi = -(a/2) to 
Tjx = + (a/ 2 ). The two half strips lie on the same line from T)i = -(z/2) 
to -(a/2) and from a /2 to z/2. Similarly, the remainder of the solid 
boundary in the 61 plane lies on the line 5 i =h - z - i(z/2) from 
Cl = 0 to Cl = “ and on the line Ci = h - z + i(z/2) from Ci = 0 
to C 1 = “ • 


To satisfy the boundary condition at the strip, a term w|j^ = w(h-X) 
is added to the potential in the transformed plane. 

The remaining boundary-value problem is solved with the aid of two 
further transformations. 

The transformation 

§2 = slnh li^ (A7) 

transforms the region under consideration in the li plane to the entire 
right half of the I2 plane, as can be seen by following the procedirre 
outlined for the first transformation. It is found that in the I2 
plane, the slot lies along the imaginary axis from -i sin (na/ 2 Z) to 
+1 sin (ita/ 2 Z) and the solid boundary including the two half strips lies 
along the imaginary axis outside of ±1 sin (jta/ 2 Z). 

The transformation 


I 


§2 +J/+ Bin® (if) 


(AB) 
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transforms the region under consideration to “the positive real half of 
the I plane, excluding the circle of radius sin (jta/2i) centered at 
the origin. The slot lies along the half of this circle in the right 
half'plane, and the complete solid boundary lies along the imaginary 
axis outside the circle. Thus, a source at the origin in the S plane 
will satisfy the condition of no flow through the solid boundary and no 
flow across the slot. Hence, the desired potential is 

<Ii = AZn (I) + wli (A9) 

Substituting equations (a 6), (A7), and (a 8) in (A9) yields 

® = AZn jsinh | (h-x) + ysinh^ ^ (h-x) + sin^ | + w(h-X) 


The constant A is evaluated by the last of equations (A5), 0 
being the equivalent of that eqiiation. 


X -00 


As X — > the potential simplifies to 


= A2n <2 sinh 


(h-x) 


X — > 


+ w (h-X) 


= AZn 


+ f (h-x)- 


+ w(h-X) 


= A y (h-x) + w (h-x) 

L 


The neglected terms are of order smaller than l/h-X. Thus, <1> will be 
zero at infinity if A = -(z/it)w so that 
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It was assumed that the velocity potential Is constant over the slot. 
The value of this constant is 


and 



at the slot 




sin 



= u at the slot 


Substituting this 


in the first of equations 


(A2) yields 


u + 



I 

It 


Zn 


sin 



at the slot 


Since it was assumed that u and w do not vary appreciably from 
slot to strip^ this equation applies everywhere at the wall and yields 


K 

Sx 


= 0 

Bxdn 


K 


for a plane wall. 



> 


j 


at the wall 


(All) 


In considering a ctirved cylindrical wall, it appears that the above 
results are not altered appreciably if the radius of curvatiire of the 
wall is everywhere large canpared to the slot spacing.^ Hence, it can be 
assumed that equations (All) are applicable to any slotted wall. 


^For a circular cylindrical slotted wall, solution of the boundary value 
problem 

59 

— = 0 at slots 

58 


— = -w at strips 
5r 


9=0 at r = 0 


yields a value for K Identical with that obtained from equation (All) . 
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In an attempt to take viscous effects in the slots into account, it 
can he assumed that as in the case of the "ideal" porous wall there is a 
pressure drop through the wall which is proportional to the normal 

velocity at the wall. In that case the first of equations (A2) is . 

replaced by (u + u) +(i/r)w = 0 at the slots which leads to 

^ + K ^ ^ = 0 everywhere at the wall . (A12) 

ox oxon R on — 

where K remains the same as in equation (All) and R is to be deter- 
mined experimentally. 

Equation (AlO) can be used to calculate- the neglected variations 
in flow quantities near the wall and also at the model if the variations 
are not negligible there. 


DISCUSSION OF A CRITICAL ASSUMPTION IN THE DERIVATION 
OF THE BOUNDARY EQUATION 

It is assTmed that the pei*turbation pressure at the wall is pro- 
portional to u + u. This is a good approximation only if 

(v + v)^ + (w + w)^«2 (u + u) U 

in addition to the usual requirements for linearization. This addi- 
tional requirement can be reasonably relaxed to 

(w + w)^ «2uU . at the walL (A13) 

Equation (AlO) Indicates singularities in w at the edges of the 
slots. Experience with wing leading edges indicates that this dis- 
crepancy can be reasonably Ignored. Howevei^_ equation (A13) should at 
least be satisfied at the center of the slots. The effect of this 
condition on u and v will now be determined. Differentiating 
equation (AlO) yields 
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Then 



and 


(w + w) 




BO that equation (AI3) hecomes 


■w^ 


« 


sin^ 



2 Uu 


or 


■w^ « 2 sln^ f — ^ Uu at the wall (Al^t-) 

\2lJ 

This result places a lower limit on the ratio of open to total area 
(a/l) for which the results of this analysis can be expected to apply 
to slotted sections. 
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APPEHDIX B 

BLOCKAGE IN A CIRCULAR TUNNEL WITH 
NONLINEAR VISCOUS EFFECT 


If, in the escperlmental determination of the porosity parameter R, 
it is found that an additional term of the form f (Scp/5n) is needed in 
the boundary equation, the result is 


^ + K + i 22 + f = 0 

Sx Sx3n R Sn 


at the vail 


(Bl) 


or 


is!! + K + iiS!! = _iSi_K 1 iSi _ f + is!!\ (-ro^ 

hx SxSn R Sn ^ ^x^n R ^n \ 3n hnj 


This type of equation cannot he solved exactly hy any presently known 
method because of the nonlinearity. However, the equation becomes linear 
if the cp^ in the nonlinear term is neglected under the assumption that 


f I 

\ Bn Bn / 



at the wall (B3) 


Hence, the equation 


Bcp* ^ ^ B^9* ^ 1 B9* 
Bx BxBn R Bn 


B9i _ B^9i 
Bx BxBn 



will be discussed. 


The transforms G* and Gi of 9* an.*! 9i remain the same as before, 
with the exception that the factor A(g) must be evaluated from the 
transform of equation (b 4) . For this purpose, the transform of 
f(B9i/Br) is needed and can be found by evaluating 



(B5) 


r=ro 
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The transform of equation (Bit-) is found to he 

A(g) - Iglo (Pgro) - igKPgli + i Pgli 1 = 

L R J 


^ ig ^-igKo + ig p [g {KKi - i p |g |Ki^ - Gf 


Upon solving for A(g) and substituting in the inverse transform 
expression for <P* equation (3^) is obtained plus the additional term 


w 


^ Pgli(Prog) - Ig lo(Prog) + KPgli(prog) 


Gflo(Prg) 


With the substitution q = Prog this becomes 


5 <p 


00 




io(<i) + ^ 1(1)' 


i _32S 

G^e" 


I 


I„(^) H (b 6 ) 


As an example of the use of equation (b 6) let 

3 


f ^ Y 

\9n/ \Sn/ 


(b7) 


This function is chosen rather than one proportional to (S*P/^n)^ 
because the pressure drop should be an odd function of the normal 
velocity. 


Hence, 


Gf = 


_ _c_ / 

' u® \ 


(3P^ro)' 


w 

if 


[ x^ + P^ro^] 


15/2 




From reference 4 it is found that 
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00 

-I ^ 


(x^ + 3®ro^) 


15/2 


igx , lei Ky (Pro]g|) 

e ox = 


( a) 

Differentiating three times with respect to g yields 


1 r -lx® _igx , 


^ r (2?rcf S« 


^ 'g| KV(P=^o|e|) 


SO that 


G-45 = - i 


U® Jl r 2-= 




or 


Of --IP ^ 2y7^j,/w^ ['^1 k,(|o|)](b8) 

Suhstituting this in equation (b 6) and separating real and 
imaginary parts yields 


5cp* = - 


(l 


mr 


w 

/■ 


®r 


■ 3 — 5 

-yfi ^ 




Iq (a) + — qli(q) 

ro 


pro 

q dq3 


=== (e^) 


f 


-,2 


] 


+ j Io(q) + iL qli(q) 


w 

1/ 



q"^ (q) 

"" VW 

, /4r\ 


;iif (.); 

2 

+ 

lo (<i) + ^ all (a) J 


(B9) 
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Upon differentiation with respect to x and setting x=r= 0 , there 
results 


6 Au = 


C / mr 3® 3 r 

d° 

Bi(q) 

q’^K^(q) 

dq 

iF Wp'^ro®/ ^ 

^ Ii(q) 

Lr J 

2 

+ 

Io(<l) + ^ Ix(q)T 
^0 


From reference 6 it is found that 


and 


qjKyiq) = (1152 + I68q2 - 5 q'‘)q^Kr(q) + ( 5?6 + 12q^ - q‘*)q®Ko(q) 


(?) 


1 2 
5005 n 


so that 
5 ^u = / 


c P 1 2 

r(ll52 + l68q2- 

5 q*)q^Ii(q)Ki(a) dq 

UP R 5005 jt ' 

0 

rs 

+ 

do(‘i) + — ^i(q) 

2 


^ |_R 


^o 



(576 + 12 q 

2 

q^) q°ii(q)K;^(q) 

[i 

2 

+ 

Io(q) + ^ Ix (q)' 
^0 

2 


dq 


(BIO) 


Setting K/ro equal to zero in (BIO) and recombining it with 
equation (37) yields 



Au — 


To 


Io('l)Ko(q) Ki(q)li(q) 



Io(<l) 

2. 

+ 

- iiCq)T 

L R J 



LO 


/ _c_ P _0^ r 

\kRTo^ J p3 E 5005 J 


■[(1152 + I68q= - 5q;*>Ii(q)Ki(q) + (5T6 +12q‘ 

• 

-q*)aii(q)Ko(q.) 3 


io(q.)l + ii(q.) 

|2 

J 


(Bll) 


aa the "blockage in a porous tunnel for vhlch the boundary equation la 

3 


^ + i ^ ^ f =0 at the wall 

ax R ^ ^ vary 


The linear-vlscouB-effect term la zero when equation (Bll) la evaluated numerically at 
P/r 1.21, and 


Au 


Mr 


2n^p^ro® 


Uo.m 

■oV fP 


courpared with the value of 


= 




Mj* 


2 rt^P®rQ® 


[2.1^1 


for a solid wall. 
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Since the factor r-^/ksccQ^ is very small, the nonlinear contribu- 
tion to the blockage correction should be negligible at subsonic speeds 
where P is not small. But as contrasted to the linear viscous effect 
which becomes small as M approaches 1 , this nonlinear viscous effect 
may become large. 
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APPENDIX C 

DERIVATION OF THE BOUNDARY EQUATION 
FOR A TAPERED SLOTTED WALL 


In this appendix an approximate smoothed boundary equation for a 
slotted wall with tapered slots will be derived. 


The development in this case is identical with that in Appendix A 
up to the point where it is found that 


cp = w 




at the slots 


or 

cp = w K at the slot 


(Cl) 


When the slots are tapered, the slot parameter K is a slowly 
varying function of x which can be expanded in a power series about 
X = 0 so that neglecting the higher order terms 


K(x) = K 



X 


X 0 


Then 


9 = K 


w + 


X = 0 


§K 

dx 


xw at the slots 


X = 0 


and 


u = ^ = K 

5x 


^ ^ dK 


X = 0 dx 


w + M 

X = 0 dx 


— X at the slot 

X = 0 


AssTjming that the last term will be negligible compared to the others, 
we have 


u = K 


w 


X = 0 


^ dK 
^ dx 


X = 0 


at the slots 


(C2) 
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Substituting this equation in equation (A2) of Appendix A yields 


u + K 


^ 

X = 0 dx 


V = 0 at the slots 


X = 0 


Then, as in Appendix A, since u and w do not vary appreciably from 
slot to strip 


Sx 


+ K 


dK 


X = 0 ^xSn dx 


= 0 Sn 


= 0 everywhere at the wall 


X = 0 


> (C3) 


K(x) = - In 


ita(x) 

2K5Ey 


sin 


Hence, the parameter 1/E in equation (6) of the text can be interpreted 
as the quantity dJC/dx | ^ which is related to the taper of the slots 

in a slotted wall having purely potential flow in the slots. 
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(a) Two-' dimensional -flow siotted test section. 



Figure A - Cross-sectional diagrams of the two slotted test sections in 
a plane perpendicular to the free stream. 




Figure 2. - Simultaneous values of slot parameter and porosity 
parameter for zero blockage in a two-dimensional-flow tunnel. 



Figure 3r Blockage factor as a function of slot parameter in a 
two- dimensional-flow, ideal, slotted tunnel. 
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Figure 4.-^ Simuitoneous values of slot parameter and porosity 
parameter for zero blockage in a circular cylindrical tunnel. 



Figure 5.- Blockage factor as a function of slot parameter in 
a circular cylindrical t ideals slotted tunnel. 




Transformed plane, ^2 


Transformed plane, S 


Figure 6.- Sketch of complex plones (Appendix A). 
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